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Abstract. We have carried out measurements of domain wall dynamics in a 
Pt/Go/GdOa,(t) wedge sample with perpendicular magnetic anisotropy. When driven 
by an easy-axis field Hz in the presence of an in-plane field Hx, the domain wall 
expansion along ±x is anisotropic, as expected for samples presenting Dzyaloshinskii- 
Moriya interaction. In the creep regime, the sign and the value of the domain wall 
velocity asymmetry changes along the wedge. We show that in our samples the domain 
wall speed vs. Hx curves in the creep regime cannot be explained simply in terms of the 
variation of the domain wall energy with Hx^ as suggested by previous works. Eor this 
reason the strength and the sign of the Dzyaloshinskii-Moriya interaction (DMI) cannot 
be extracted from these measurements. To obtain reliable information on the DMI 
strength using magnetic field-induced domain wall dynamics, measurements have been 
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performed with high fields, bringing the DW close to the flow regime of propagation. In 
this case we find large values of DMI, coherent with those obtained from current-driven 
domain wall dynamics. 


PACS numbers; 75.70.Ak, 75.60.Ch, 75.60.Jk 


Chiral magnetic textures such as Dzyaloshinskii domain walls (DDW) [1] and 
skyrmions [2] are attracting attention because of their possible applications as 
information carriers in spintronics devices. DDW are Neel walls with a fixed chirality, 
stabilised, in non-centrosymmetric stacks, by the Dzyaloshinskii-Moriya interaction 
(DMI) [3, 4] present at the interface between a magnetic layer and a heavy metal with 
large spin-orbit coupling. When driven by a Spin Hall effect related spin-orbit torque 
(SHE-SOT) [5, 6, 7] DDW in systems with perpendicular magnetic anisotropy (PMA) 
move with large efficiency [8, 9, 10]. Afso, it has been predicted that isofated skyrmions 
injected in nanotracks can be moved with very low current density and are moreover 
insensitive to defects [11]. Engineering materials with large DMI has therefore become 
an important issue both for domain wall and skyrmion physics. 

So far ah-initio calculations of interfacial DMI are rare and concern perfect 
interfaces difficult to compare with the mixed interfaces found in “real” samples 
[12, ?]. The information presently available on the DMI strengths relies on experimental 
work. A large input has been given by Spin-polarised Scanning Tunneling Microscopy 
measurements that show the presence of chiral magnetic textures or skyrmions in 
systems consisting of one monolayer of Fe (or Mn) on heavy metal substrates 
[13, 14, 15, 16] in ultra-high vacuum and at low temperature. In the last few years, 
domain wall dynamics and nucleation measurements at room temperature have revealed 
the presence of DMI in less ordered, non centrosymmetric ultrathin magnetic layers 
with PMA, made by magnetron sputtering [9, 10, 6, 17]. More recently, Brillouin 
light scattering experiments have also highlighted the presence of DMI in similar PMA 
samples [18, 19]. 

It has been shown recently that when, in a nanostrip or in a bubble domain, an easy- 
axis field drives the DW dynamics in the presence of an in-plane field (aligned 
along +x), the DW speed is different for up/down and down/up DDWs propagating 
afong ±a: [20, 21, 5]. This phenomenon is related to the symmetry breaking introduced 
by the in-plane field. The Dzyaloshinskii-Moriya interaction acts as a longitudinal 
chiral field T^dmi = D/ {iiqMsA) (where D is the DMI strength. Mg is the saturation 
magnetisation and A is the domain wall width parameter) localised on the domain walls, 
having opposite directions for up/down and down/up DWs. Beyond a critical strength, 
the DMI forces the DW magnetisation in the Neel configuration (see sketch in Figure 1) 
[1]. Although the in-plane field does not drive the dynamics, it will respectively stabilise 
{vs. destabilise) the DWs having their magnetisation m parallel (vs. antiparallel) to it. 
For a parallel {vs. antiparallel) alignment between and m the DW speed increases 
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{vs. decreases) with respect to the Hx=0 case. In the high speed (flow) regime, the speed 
increase (vs. decrease) is mainly due to the widening (vs. narrowing) of the DW with 
[5]. In the low speed (thermally-activated or creep) regime, the speed dependence 
on Hx has been related to the variation of domain wall energy [20]. The DDW width 
(resp. DW energy) is expected to have a minimum (resp. maximum) value when the 
applied in-plane field is equal and opposite to the stabilising i^oMi field i.e. when the 
DW acquires a Bloch form. In the two DW propagation regimes, this is the field 
for which the DW speed is predicted to exhibit a minimum. With these assumptions, 
Hx could therefore be a direct measure of the DMI energy density D, provided that the 
domain wall width parameter A = sjA/K q (Kq being the effective uniaxial anisotropy 
and A the exchange constant) and Mg are known. 

In the following, we will show that the DW speed vs. in-plane held curves in the 
creep regime cannot in general be used to extract the strength and the sign of the 
DMI, as was done for Pt/Co/Pt samples [20, 21]. Moreover, we find that the v{Hx) 
curves measured for the same sample in the thermally activated and in the flow regimes 
can have different trends. Although the mechanism determining the exact trend of 
the velocity curves in the creep regime is not clear, we show that it cannot always be 
described simply in terms of the variation of DW energy with H^. Our measurements 
on Pt/Co/GdOa; films suggest that modifications of the pinning barrier landscape upon 
application of the in-plane field also contribute to the trend of the v{Hx) curves. 

A Pt(5nm)/Co(lnm)/Gd(t) stack with varying Gd thickness {t = 2 — 5nm) was 
grown on a Si/Si 02 substrate by magnetron sputtering in the shape of a wedge, and 
oxidised by O 2 plasma for 35 seconds. Gonsequently 2nm of A1 were deposited on top of 
the stack to protect it from further oxidation. The varying thickness of the Gd layer is 
at the origin of a gradient in the oxygen content at the Co/Gd interface, which varies the 
interfacial anisotropy [23]. All the samples present a well defined PMA, with in-plane 
saturation fields varying between 1.6T (for 2nm Gd) and 0.6T (for 5nm Gd). Domain 
wall dynamics was studied at room temperature by wide-field Magneto-Optical Kerr 
microscopy, using a combination of easy-axis and in-plane magnetic fields. Hz pulses 
of amplitude ~10mT and duration ~20-100ms were obtained using a conventional, 
uncooled coil. The Hz pulses, driving the displacement of the DWs, were applied in 
the presence of a continuous in-plane field H^, along ±x, which tunes the stability of 
the DDW internal structure. With such amplitudes of the Hz field, DW speeds are of 
the order of some O.lmm/s, the dynamics is thermally activated and described by the 
so-called creep regime. 

Starting from (down or up) saturation, a bubble domain was created by applying 
an up or a down Hz pulse. The image of the domain was saved as a reference image. 
An Hz pulse was then applied to enlarge the domain by DW propagation, and the 
new image was acquired. The difference between the two images gives the domain wall 
displacement that occurred during the field pulse. A black (white) contrast in the images 
corresponds to the expansion of an up (down) domain. The domain wall speed in a given 
direction can then be extracted from the ratio of the DW displacement and the pulse 
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Figure 1. Left; Expansion by DW propagation of an up (black contrast) and a down 
(white contrast) domain in samples (A), (B), (C), (D). The Gd thickness increases 
from 2nm to 5nm going from (A) to (D). The DW displacements are obtained by 
an Hz field pulse with amplitude ~10mT and duration ~20-100ms and a continuous 
in-plane field = -1-200 mT. Right; schematic view of a bubble domain expansion; 
the red arrows represent the equilibrium orientation of the magnetisation at the center 
of the DWs. 


duration. DW displacements in the ±.x directions were measured for a fixed value of 
the Hz field, for several values of the in-plane field between -300 mT and +300 mT. In 
order to correct the residual Hz component that may arise from a misalignment of the 
in-plane electromagnet, measurements were taken for both down and up domains. 

Figure 1 shows the differential images recorded in four positions of the wedge sample 
(called samples (A) to (D) from now on) corresponding to increasing values of the Gd 
thickness (from 2 to 5nm) for Hz field pulses of the order of 10 mT and an in-plane 
field of +200 mT. Without in-plane field, the propagation of the DWs is isotropic and 
the domains are round. Similar to previously reported experiments, the H^ field breaks 
the rotational symmetry and the propagation becomes asymmetric in the ±a: directions. 
Note that the sign and the amplitude of the speed asymmetry depend on the sample 
composition. Indeed, in sample (A) the down/up DWs move faster than the up/down 
DWs while in sample (B) the asymmetry is practically vanishing, i.e. up/down and 
down/up DWs move at the same speed. In sample (C) the DW speed asymmetry 
reverses with respect to (A), i.e. the up/down DWs move faster. Finally, in (D) the 
asymmetry found in (A) is recovered. 

According to previous work [20, 21], the cancellation (resp. change of sign) of the 
DW speed asymmetry may be attributed to a vanishing (resp. reversed) value of the 
DMI. This result is unexpected and counter-intuitive. As one moves across the sample, 
from (A) to (D), the decreasing degree of oxygen content modifies the composition of the 
Co/Gd interface, as shown experimentally by the changing PMA. However the sample 
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presents a considerable PMA for the thinner Gd layers, which is an indication that the 
oxidation only concerns the top Co interface. Therefore the bottom Pt/Co interface, 
which is expected to provide the most important contribntion to the DMI [12], should 
not be affected by the varying Gd thickness. This is confirmed by X-ray reflectivity 
data. 

In order to clarify the interpretation of the DW dynamics in the creep regime and 
to have an independent measurement of the sign of the DMI, we carried out current- 
induced DW dynamics measurements. For this purpose, the samples were patterned 
into 1 pm wide strips by e-beam lithography and the DW dynamics was studied for 
a fixed value of the current-density J = 1.2 x lO^^A/m^ and variable values of Hx- 
The results show that for all samples (note that sample (C) could not be measured, 
due to deterioration during the patterning process) the domain walls move in the same 
direction, opposite to the electron flow. Since in these systems the direction of the DW 
displacement is determined by the sign of the Spin-Hall angle in Pt (which is the same 
for samples (A) to (D)) and by the chirality of the DDW [1, 9, 10], this results is a 
strong indication that the domain walls in all the samples have the same chirality and 
therefore the sign of the DMI is sample independent. The results of the current-driven 
DW speed vs. Hx field curves for samples (A) and (B) are shown in Figure 2 . 



Figure 2. Left: Domain wall speed vs. curves measured with constant current 
density J = 1.2 x 10^^ A/m^ for samples (A) and (B), for which a large and a vanishing 
field-induced domain wall speed asymmetry are found in the creep regime respectively. 
Right: Differential Kerr image showing an example of the displacement of DWs in 
nanostrips. 


The speed variation as a function of in-plane field is similar to that shown by 
other authors in strips of DMI materials [9, 10, 24]. In all the curves, the speed of the 
down/np DWs increases for positive Hx fields and decreases for negative fields. The 
symmetric curve is found for the np/down domain walls, as expected for chiral Neel 
walls. If we neglect the rotation of the magnetisation within the domains, the domain 
wall speed driven by the current J via the SHE-SOT can be expressed as [1] : 

A TT 


^ = 7o —-xM^cos^ 

a 2 


( 1 ) 


where 70 is the gyromagnetic ratio, a is the damping parameter, A is the domain 
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wall width, if is the angle of the DW magnetisation with respect to the x-axis, and 
X = hOnJ/{2efioMft) where 9h is the Spin Hall angle and t the magnetic layer thickness. 
It can be shown (see Suppl. Information) that for our samples the variation of cos'ip 
with is negligible except around = —i^oMi where it changes sign, so that the 
v(Hx) shape is mainly determined by the modification of the domain wall width with 
Hx- Since the DW width increases for an field parallel to the DW magnetisation, 
our measurement show that down/up DWs have their magnetisation parallel to the +x 
direction and therefore that the DWs in the Pt/Co/GdOa, samples have left-handed 
chirality, like in Pt/Co/AlOa, [17, 5]. This is not surprising, as we expect that the DMI 
interaction is mainly located at the Pt/Co interface. 

The velocity of the down/up DW in sample (A) changes direction under the effect 
of a negative in-plane field ~ —280mT; this is associated with the switching of the 
DW chirality when the negative field exceeds the local chiral T^dmi field. This in¬ 
plane field value is therefore a measure of T/dmi- Note that in sample (B) the switching 
of the DW velocity is hindered by the larger DW pinning [9, 24]. 

The constant sign of the DMI for all the samples - assessed by the constant direction 
of current-driven DW motion at zero field - in contrast with the different DW 
velocity asymmetries observed for the different samples in Figure 1, sheds doubts on 
the possibility to deduce the sign of the DMI from the domain expansion images in the 
creep regime. In order to clarify the interpretation of the field-induced measurements, 
we measured the DW speeds as a function of field for the bubble domains shown in 
Figure 1. 

The velocity curves are shown in Figure 3 for the two domain walls propagating 
along the x-axis and having their magnetisation either parallel or antiparallel to the 
field. The up/down and the down/up DWs exhibit the same behaviour for opposite 
fields, as expected for chiral Neel walls. The curves for sample (D) - corresponding to 
the thicker Gd layer - present the main features found by other authors for DDWs in 
Pt/Go/Pt films [20, 21 ]. The speed of the down/up DW increases for a positive in-plane 
field, and for negative fields it decreases down to a minimum value between -100 mT 
and -200 mT, where the velocity starts increasing again. On the other hand, the curves 
measured for samples (A-G) strongly deviate from the expected behaviour, showing in 
particular a maximum rather than a minimum in the DW speed. 

In the thermally activated regime, the DW velocity is given by: 

v(H^) = Voexp{-r]Hf>^) ( 2 ) 

where vq is the characteristic speed, /r=l/4 is the creep scaling exponent and rj = 
UcdiV^/ksT where Uc is an energy scaling constant and i/cut the critical magnetic field 
[25, 26]. Following Ref. [26], Uc is related to ^ (the correlation length of the pinning 
potential) and to the Larkin length Lc = characteristic length of 

rigid microscopic DW segments) and iLcUt = where udw N the DW energy 

and 7 is the pinning strength of the disorder. By assuming that neither ^ nor 7 are 
modified by H^, Je et al. [20] conclude that the shape of v(Hx) is solely due to the 
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in-plane field dependence of the DW energy. According to [17], the energy of a DDW, 
taking into account the modification of the DW profile with reads: 

_ 2 D 

a = (Too[Vl - h^ + {h + )(arcsinh =F 7r/2)]. (3) 

TT L>cO 

where ctoo = 4:^/AKq is the DW energy at rest, Dco = 4^1 AKq/tt = ctoo/tt gives 
the onset of magnetisation cycloids, h = Hx/Hko and the =F signs refer to the DW 
having its magnetisation parallel/antiparallel to the field. The energies of the DW 
favoured/unfavoured by the in-plane field are the same when h = —(2/7r)(D/Dco) or 
when Hx = —D/(fj,QMsA) = -Hbmi- This is the in-plane field for which the DW energy 
is maximum. From equation 2 it then follows that the DW velocity should exhibit a 
minimum for Hx = —T/dmi- This is indeed observed for sample (D). Note that the 
left-handed DW chirality deduced from the measurement agrees with the results of the 
current-induced measurements. 
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Figure 3. Domain wall speed vs. field measured in the thermally activated regime 
for bubble domains in Pt/Co/GdOa^ samples (A) to (D), for the DW propagating along 
the x-axis direction. 


The speed vs. Hx curves obtained for samples (A-C) show a different behaviour. For 
sample (A) the speed asymmetry is the same as for sample (D), but the velocity of the 
down/up DW increases for negative fields and decreases up to a critical field for positive 
fields. In sample (B) the speed asymmetry practically disappears and the speeds of the 
up/down and down/up DWs continuously decrease with both positive and negative Hx 
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fields. In sample (C) the asymmetry is switched for down/up and np/down DW, with 
respect to sample (D). Therefore in these three samples the v{Hx) curves do not follow 
the variation of the DW energy. Curves deviating from the expected behaviour have 
also been recently reported in the literature [27]. 

Note that the “anomalous” curves are found in particular for samples (B) and (C), 
for which the sign of the speed asymmetry would suggest that the value of D is either 
vanishing (for (B)) or opposite (for (C)) to the one of sample (D). This indicates that 
in the creep regime extreme care should be taken when extracting information on the 
DMI sign and amplitude simply on the basis of the asymmetry (or lack of asymmetry) 
of the Kerr microscopy differential images. Before assessing about D, the full speed vs. 
Hx curves should be examined and compared with the curves predicted by the existing 
theoretical models. 

In order to verify the role of the DW pinning on the speed vs. field, we have 
repeated the field-dependent measurements for larger values of the Hz fields, bringing 
the domain wall velocities to a regime (^ Im/s) where the propagation is much less 
sensitive to the pinning generated by local variations of the anisotropy field. Pulsed 
Hz fields up to 200 mT and duration down to 20 ns were obtained using a 50 p.m wide 
microcoil coupled to a fast current pulse generator [28]. The results reported in Figure 4 
for samples (A) to (D) show that in these conditions the speed vs. H^ curves all acquire 
the trend expected for chiral Neel walls in the flow regime [5]. 

In the high field regime, the stationary DW velocity is given by n = 'jQATHz/amzo, 
where At is the Thiele domain wall width [2] and uizo is the easy-axis magnetisation 
within the domains. The speed variation with H^ is mainly related to the modification of 
the Thiele DW width with the in-plane field (see Suppl. Information). In all the samples, 
the down/up DWs propagate faster than the up/down DWs for positive H^ fields, 
confirming once again that the DW chirality is the same (left handed) in agreement 
with the current-induced measurements and the field-induced (creep) measurements for 
sample (D). For samples (A) to (C) the DW speed of the down/up DWs decrease down 
to the largest available negative H^ field, with a saturation but not a clear minimum in 
the DW speed. This suggests that T^dmi in these samples is of the order or more than 
+300 mT. 

For sample (D), where the PMA (and therefore the T^dmi field) is reduced, the 
down/up DWs exhibit minimum speed for H^ =~ —180mT, a value close to that found 
for the same sample in the creep regime. By taking Mg = 1 x 10® A/m (measured 
by VSM-SQUID), HqHk = 0.7T (measured by EHE) and A = 2.2 x 10“^^ J/m [30] 
the expression T^dmi = D/(fj,oA Mg) gives rise to a value oi D = 1.27mJ/m^, with 
A = 7.1 nm. Taking into account the larger thickness of the Co layer (Inm) in 
our Pt/Co/GdOa; samples, this value scales reasonably well with the D = 2mJ/m^ 
value found for Pt/Co(0.6nm)/AlOa, [17, 5]. For samples (A) to (C) it is difficult to 
obtain a precise value of D from the field-dependent measurements, where the minimum 
speed is not well defined. The value of D for sample (A) may be derived from the in¬ 
plane field for which the DW chirality switches when driven by spin-polarised current 
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—280 mT). Using the values of Mg and A used for sample (D), and the measured 
in plane saturation field iiqHk = 1.6 T giving rise to A = 4.7nm, we obtain a value 
oi D = 1.31 mJ/m^. We estimate that the uncertainty associated to the value of the 
exchange parameter A, together with the error associated to the definition of the 
field where the DW velocity is minimum, allow the determination of D with a precision 
not better than ±0.2mJ/m^. The similar D values found for the two samples indicate 
that the DMI strength is homogeneous along the wedge sample and that it is mainly 
arising from the Pt/Co interface. 

As a consequence, the “anomalous” v{Hx) curves in the creep regime do not bear 
any information about the sign and strength of the DMI. In Figure 3, the value for which 
the speed is maximum in samples (A) and (C) is not related to the D value, and the 
absence of speed asymmetry for sample (B) is not a signature of a vanishing D. Since 
the anomalous behaviour of the v{H^) curves is observed only in the creep regime and 
for samples (A) to (C), we conclude that this feature may be related to modifications of 
the domain wall pinning with which depends on the details of the Co/Gd interface. 



Bx (mT) Bx (mT) 


Figure 4. Domain wall speed vs. field measured for bubble domains in 

Pt/Co/GdOa; for samples (A) to (D) for 20ns-long Hz field pulses varying between 
70mT and 200mT. For samples (B) to (D) the curves were measured for two Hz 
field values (empty symbols correspond to the scale to the right). The trends of the 
normalised speed curves are the same for each field value. For sample (D) the speeds 
are larger, as the depinning of the DWs occurs for lower fields. 
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Since the measurements were taken with ms-long pulses in the creep regime and with 
ns-long pulses in the flow regime, the effect of the pulse length on the DW pinning may 
also play a role. 

Some information on the nature of the top interface, as the presence or not of 
CoO, can be obtained from the temperature dependence of magnetic hysteresis loops. 
We have carried out magnetisation measurements with variable temperature between 
10 K and 300 K in a VSM-SQUID of Quantum Design (Figure 5). For sample (C) a 
change of the hysteresis loops, which are square with 100% remanence at 300 K, is 
observed around 225 K, where they become partly tilted and the remanence decreases 
to about 60%. This indicates a decrease of the PMA. Upon decreasing the temperature 
further, the coercivity increases strongly and below 70 K a shift of the hysteresis loop 
to negative fields develops. Both observations can be attributed to the presence of an 
ultrathin layer of CoO at the Co/Gd interface, which becomes antiferromagnetic around 
225 K with a blocking temperature around 70 K. For sample (D), the only one presenting 
“expected” v{Hx) curves in the creep regime, the cycles do not exhibit any exchange 
bias indicating that no CoO is formed at the top Co interface. 
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Figure 5. VSM-SQUID measurements carried out from 10 K to 300 K for an out- 
of-plane field up to 1 T. Left: in sample (C), the shift of the cycle at low temperature 
is an indication of the presence of CoO at the top Co/Gd interface. Right: in sample 
(D), the cycle does not exhibit a shift, sign of the absence of relevant oxidation. Note 
the factor 10 difference in the field scale. 

The “anomalous” behaviour of the v{Hx) curves in the creep regime seems therefore 
to be related to the presence of Co oxide at the top Co interface, and the details of the 
curves to the different degree of oxidation. Although the CoO is paramagnetic at the 
room temperature, it exhibits a magnetic susceptibility in the x-direction [31]. We 
speculate that the CoO magnetic moments induced in the x direction by the in-plane 
field may act as an extra pinning potential acting on the DWs. Since the magnetic 
susceptibility may depend on the CoO thickness, this could explain why different samples 
exhibit maximum velocity for different fields. The description of the creep law simply 
in terms of the variation of the DW energy may not be general, as the pinning potential 
landscape may be strongly affected by the in-plane field. 

In conclusion, we have shown that in Pt/Co/GdO^ samples with different oxidation 
degrees of the Co/Gd interface the dependence of the DW velocity as a function of 
the in-plane field cannot be interpreted within the creep law relating the DW speed 
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changes exclusively to the DW energy variations. Therefore in these samples, the v{Hx) 
curves fail to give information about the sign and the strength of the DM interaction. 
We have correlated the failure of the proposed creep law with the modification of the 
pinning potential landscape induced by the in-plane field. When by applying strong 
and ultrashort ont-of-plane field pulses we change the dynamic regime of the DW 
propagation, the v{Hx) curves indicate that the chirality of the DDW is left-handed, 
and D is of the order of 1.3mJ/m^ for Inm Co. 

This work was supported by the Agence Nationale de la Recherche, project ANR 
11 BSIO 008 ESPERADO. SP acknowledges the support of E. Wagner and of the staff 
of the Nanofab facility in Inst it ut Neel. 
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SUPPLEMENTAL MATERIAL 


1. MICROMAGNETIC SIMULATIONS 


Micromagnetic simulations of domain wall (DW) velocities driven either by an out-of- 
plane magnetic field or by an electrical current in a nanostrip were performed using a 
full ID micromagnetic model, as introduced in [1], The unknown parameter is the full 
m(T, f) profile. The demagnetizing field is evaluated by direct summation for a given 
nanostrip width (500 nm). In our case, as the sample thickness and the DW width are 
much smaller than the typical strip width, the wall may be considered as infinitely long. 
The computed value of the demagnetizing field is averaged over the sample thickness 
and over the strip width. In practice, only the x component of the demagnetizing field 
is evaluated; the y component that is transverse to the strip is negligibly small and the 
z component is approximated by a local value with unity demagnetizing factor. For 
the evaluation of DW dynamics, the finite box containing the DW is shifted along the 
strip so that the DW is always kept in its center. For field-driven dynamics, the case of 
bubble domains in the continuous layers is too complex to be treated specifically; the 
DW propagation along a given direction has been assimilated to that of a domain wall 
within a strip aligned along that particular direction. Having obtained the complete 
profile of magnetization in the domains and across the DW, all quantities of interest can 
be evaluated. These contain: 

(i) the Thiele DW width [2], defined as: 


A = 1 /■ 

At S J dx 


( 4 ) 


where S is the cross-section area of the nanostrip, oriented along the x direction. 

(ii) the angle $ of the DW magnetic moment which controls the force exerted by 
the Spin Hall Effect (SHE) on the DW in the case of current-driven DW motion. In 
the presence of an in-plane field and SHE the domain magnetization rotates so that the 
straightforward evaluation of cos $ fails. To remedy this, we evaluate by integration the 
SHE force on the DW. This is proportional to f (m x dxrh)y, with a value of vrcosd* 
in the simple case, tt standing for the angle between the two domain magnetizations. 
Thus, to compute d>, this integral is numerically evaluated, and then divided by the 
angle between the domain magnetizations. 

The micromagnetic parameters chosen for the simulations are Mg = 1000 kA/m, 
A = 22 pJ/m, Ku = 0.87 and 1.17 MJ/m^ where Kq = Ku — iJ,oMf/2 is the 
effective anisotropy including the perpendicular demagnetizing field effect in the local 
approximation. Eor the strength of the Dzyaloshinskii-Moriya interaction we used 
the value D = 1.8 mJ/m^, which is slightly larger than the one evaluated for our 
samples. Eor the magnetization dynamics, the gyromagnetic ratio of the free electron 
7o = 2.21 X 10^ m/(A-s) and the damping factor a = 0.5 extracted from DW dynamics 
experiments were used. 
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The aim of these simulations is to explain qualitatively the observed trends in the 
speed vs. curves, in particular for samples with different magnetic anisotropy values. 

1.1. Current-driven dynamics 

The simulations were carried out for a fixed value of the current density J = 1 x 10^^ 
A/m^. Figure 1.1 (a) shows the variation of the domain wall speed as a function of 
longitudinal in-plane field H^, for two values of the out-of-plane-anisotropy Ku = 0.87 
and 1.17 MJ/m^ (A and B). The figure shows that the DW speed changes more rapidly 
for the low anisotropy value, in agreement with the experiments. For strong enough 
negative field (i.e. antiparallel to the DW magnetisation direction) the domain wall 
velocity changes sign, due to the reversal of the domain wall magnetisation. The reversal 
field (Hx = —Hdmi = —D/(/ioAMs)) is smaller for the smaller anisotropy value, as the 
domain wall width is larger. As seen in the main text, the detailed shape of the speed vs. 
Hx curve depends both on the variation of the Thiele domain wall width and on the 
value of cos$, where d> is the angle that the DW magnetisation forms with the a:-axis 
along which the in-plane field is applied. Figures 1.1 (b-c) illustrate the simulated values 
of At and cos $. For positive values the speed curve variations are only related to 
the variation of A^, which is larger for the smaller anisotropy. For negative fields, the 
change of the cos $ sign (reversal of the DW magnetisation direction) is at the origin of 
the reversal of the DW speed. 
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Figure SI : (a) Current-driven domain wall speed vs. (IqHx calculated for a fixed 
value of the current density J = 1 x 10^^ K/w? and two values of the anisotropy energy 
Ku = 0.87 MJ/m^ (curve A) and 1.17 MJ/m^ (curve B); (b) variation of Thiele DW 
width with HqHx ; (c) variation of cos$ with jioHx. 
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1.2. Field-driven dynamics in creep and flow regimes 

The field-driven domain wall dynamics under in-plane field is dependent on the working 
regime (creep or fiow). The dynamics of chiral Neel domain walls in the presence of an 
in-plane field has been studied in the creep regime [3, 4] and in the fiow regime [5] . 
It has been shown that the curves can be explained in terms of the variation of 

the DW energy with the field. Figure 1.2(a) shows the variation of the domain wall 
energy vs. for a sample having the magnetic parameters given above and fiQHz=50 
mT . The DW energy is maximum for the = —Hjjmi, where the DW has the Bloch 
form. This field is of course dependent of the anisotropy value. 
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Figure S2 : (a) Domain wall energy vs. simulated for a fixed value of 50mT. 
(b) variation of DW domain wall speed vs. assuming a vq value of 8 x 10“^ m/s 

and rj = a(Hx)/( t{0) in equation (2) of the main text. (A) and (B) refer to the two 
anisotropy values. 

The expected dependence of the DW speed in the creep regime is shown 
in Figure 1.2(b) for the two anisotropy values. The curves present a minimum for 
Hx = —Hdmi and have a symmetric behaviour on either sides of this field. The 
behaviour measured for sample (D) is in qualitative agreement with these curves. That 
found for samples (A) to (C) is in strong disagreement, showing that the v(Hx) curves 
are not simply related to the change of the DW energy with H^. 

If we neglect the tilt of the magnetisation in the domains, in the fiow regime the 
stationary DW velocity is given by n = and the speed variation with is 

expected to be related to the modification of with the in-plane field. The two curves 
simulated for a, fi^Hx = 50 mT are shown in Figure 1.2. As in the case of current-driven 
dynamics, the Thiele DW width increases (decreases) for positive (negative) in-plane 
field, and in a larger extent for the low anisotropy value (A). The dip in the DW speed 
in the vicinity of = —Hdmi, corresponding to the precessional regime of the DW, 
is not observed in the experimental data. Although this discrepancy is still the object 
of our studies, we believe that the local changes of the DMI strength over the region 
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swept by the DW during the measurements may contribute to the smoothening of the 
discontinuity associated to the precessional regime. 



Figure S3 : (a) Domain wall speed vs. in the flow regime simulated using jiffHx 

= 50 mT and the two anisotropy values given above; (b) the variation of the Thiele 
DW width. 
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